Magnetic phases of the mixed-spin J\ — J 2 Heisenberg model on a square lattice 
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We study the zero-temperature phase diagram and the low-energy excitations of a mixed-spin 
(Si > S2) Ji — Ji Heisenberg model defined on a square lattice by using a spin-wave analysis, 
the coupled cluster method, and the Lanczos exact-diagonalization technique. As a function of the 
frustration parameter J2/J1 (> 0), the phase diagram exhibits a quantized ferrimagnetic phase, a 
canted spin phase, and a mixed-spin collinear phase. The presented results point towards a strong 
disordering effect of the frustration and quantum spin fluctuations in the vicinity of the classical 
spin-flop transition. In the extreme quantum system (Si, S2) = (1, \), we find indications of a new 
quantum spin state in the region 0.46 < J2/J1 < 0.5. 
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I. INTRODUCTION 

There has recently been an increasing interest in 
Heisenberg spin systems exhibiting quantum phases with 
an extensive magnetic moment. An intriguing example 
is the bilayer Hall system at filling factor v = 2 which 
was shown to possess a cantcd,.spm state with sponta- 
neously broken U(l) symmetry.tra A closely related phase 
diagram was studied in the framework of the bilayer 
quantum Heisenberg model subject to a perpendicular 
magnetic fieldfl Some specific properties of the quantum 
phase transitions in systems with quantum ferromagnetic 
phases have been previously discussed by using a spe- 
cial class of-, lattice models with quantum-rotor degrees 
of freedoms In the present paper we study a J± — J 2 
Heisenberg model which is a natural (mixed-spin) exteru 
tion of the well known antiferromagnetic J\ — J 2 modelu 
and exhibits most of the quantum magnetic phases found 
in the aforementioned rotor systems. 

The model is defined by the Hamiltonian 

H = J\ Si r • S2 r / + J 2 (Si r • Si r / + S2 r • S2 r ') , 

(r,r') [r,r'] 

(1) 

where (r,r') and [r,r'] denote pairs of nearest and next- 
nearest (diagonal) sites of the square lattice. We choose a 
checkerboard arrangement for the spin operators Si r and 
S 2r [Si 2 = Si (Si + 1), S 2 2 = S 2 (S 2 + 1), and Si > S 2 ] 
and introduce the parameters a = J 2 /Ji (Ji, J 2 > 0) 
and a = S1/S2. With some minor modifications, e.g., 
introducing a spatial anisotropy in the nearest-neighbor 
exchange interaction, the Hamiltonian (|l|) could describe 
a large class of real mixed-spin compounds such as the 
bimetallic molecular magnets.Q 

In the classical limit the phase diagram of the mixed- 
spin system contains the ferrimagnetic state (F phase), 
the canted state (C phase), and the mixed-spin collinear 



statoo (N phase) shown in Fig. |l|. The latter phases are 
stable, respectively, in the parameter regions a < a c \, 
a c i < a < 0.5, and a > 0.5. The classical F-C tran- 
sition at cv c i = (2cr) _1 is continuous, whereas the C-N 
transition at a = 0.5 is connected with a flop of the Si 
and S 2 spins. The canted spin phase appears as a result 
of the magnetic frustration without explicit breaking of 
the SU (2) symmetry (in the mentioned bilayer systems 
the C phase is generated by a perpendicular magnetic 
field) . Note also that the S 2 spins remain collinear in the 
C phase. The angle 8 measuring the local orientation of 
the classical Si spins in respect to the global magnetiza- 
tion axis z alternatively takes the values 



8 = ± arccos ■ 



1 

2aa 



(2) 



Bimetallic molecular magnets are normally characterized 
by small exchange constants J\ such that the canted 
phase may appear for a moderate magnetic frustration 
provided that the parameter a is large enough. In this 
article we study the quantum phase diagram of the model 
(|l|) by using the standard spin- wave theory (SWT), the 
coupled cluster methodB (CCM), and the Lanczos exact- 
diagonalization (ED) technique. 
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FIG. 1. The classical ferrimagnetic (F), canted (C), and 
collinear (N) spin phases of the mixed-spin Ji — J2 model on 
a square lattice. 
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II. PHASE DIAGRAM 

A. Spin-wave analysis 

To perform a SWT analysis we assume that the classi- 
cal spins lie in the xz plane. Using as a quantization axis 
the local orientation of the classical spins (the z' axis in 
Fig. [l]) , the leading spin- wave terms in the expansions for 
the spin operators Si r and S2 r read 



Q. 

o,, 



' {Si 



'St t 



Sil = smdi_ b f r (Si - a\ r a ir ) + cos# i: 



4r +«ir), 



■a*) (« = 1,2), 



where the angle 8i T measures the local deviations of z' 
from the magnetization axis z. A substitution of the lat- 
ter expressions in Eq. (|l|) yields the following spin-wave 
Hamiltonian written in terms of the Fourier transforms 
ciik of the boson operators <Zi r : 



i=l k 



E 

k 



, f , -^ ik I t t i \ 

Aiua^aik + -^-(a- k a-_ k + Oikai-k) 

+ E . (3) 



[Cka lk a 2-k + Acai k a2k + h. 



Note that Eq is the energy of the classical spin state; the 
coefficients in Eq. (|J) are defined below. The wave vector 
k = fciei + &2e2 runs in the paramagnetic Brillouin zone 
(PBZ) |fci ± fc 2 | < Tr/ao (the unit vectors ei _L e 2 point 
towards the main directions on the square lattice). We 
shall henceforth use a system of units where the lattice 
spacing ciq = 1. 

The canonical transformation 



a ik 



2 

E 

3=1 



(< = 1,2) 



recasts Hq to the diagonal form 



Wo = ^ -E Ak+ EE^( k ) 



i=l k 



"Ik"* + 2 



(4) 



where A k = (Aiu + A 2k )/2. In the PBZ any of the dis- 
cussed magnetic phases is characterized by two branches 
of spin- wave excitations [u>i — Wj(k), i = 1,2] 



wi,2 = V F k T vG k 



(5) 



where F k = (e? k +e2 k )/2+I>£-C£, G k = (e? k -el k ) 2 /4+ 
(e?k + 4)Pk - Ck) + 2(A lk A 2k + B lk i?2k) (Cg + £> k ) - 
4(A lk B 2k + A 2kJ B lk )C k -D k , and e lk = {A^ - BgJ 1 / 2 



(i = 1,2). and w 2 are the positive solutions of the 
biquadratic algebraic equation 



A(w) = det V(uj) = . 
The dynamical matrix T>(u) reads 



(6) 



Ai k - lu Bi k £> k 

^i k Aik + Ck 

£>k Ck A 2k - w 

Ck ^k £?2k 



C k 
S 2 k 

^2k 



A straitghforward calculation gives the following ex- 
pression for the on-site magnetizations m\ = (Si T ) and 
m 2 = (S 2 *'): 



where 



Ik 



N 4-^ 



.,(!) 



,(2) 



= 1,2), (7) 



Ai 2 (wi) 



;f)A u (wi) 



y 2k 



^20^2 - Wi)An(wi) ' 



(2) 
'4k 



(2) 
"2k 



-A n 



-w 2 ) 



2w 2 (w| - wf) ' 
-A 13 (-w 2 ) 



^-2^2(^1 - w 2 )A n (-w 2 ) 



(8) 



In the above expressions, Ay (w) is the minor of the ij 
element of A(w) and N is the number of sites of the lat- 
tice. 
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FIG. 2. The on-site magnetizations mi = (Sir ) and 
1T12 = (52^ ) plotted as a function of J2/J1: SWT (full lines), 
CCM (filled circles), ED for N — 20 spins and peroidic bound- 
ary conditions (open circles). The ED data for mi in the C 
phase are obtained from the relation mi = (Si J)/ cos # and 
the pitch angle 6 is extracted from the correlation function 
K Sl in Eq. (§). 
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B. Magnetic phases 

The spin-wave results for the on-site magnetizations 
mi and rri2, Eq. (0), in the extreme quantum system 
(1, |) are presented in Fig. |[ We suppose that in the 
ferrimagnetic phase the sublattice magnetizations are ori- 
ented along the global z axis, i.e., mi = (Si*) and 
ni2 = —(S2r)- The F phase is also characterized by the 
net ferromagnetic moment per cell Mq — (Si — S2) = 5 
(see Fig. ||). In the general case it is straightforward to 
observe that M takes only integer or half-integer val- 
ues in the F phase. Following the authors of Ref. ^, this 
state may be reffered to as a quantized unsaturated fer- 
romagnetic phase. The ED data from Figs. ^ and [| show 
that the position of the continuous phase transition a c i is 
slightly changed by quantum fluctuations towards larger 
a. a c i is fixed as a point where the spin of the ground 
state changes from (S 1 -S 2 )N/2 to (Si- 5 2 )iV/2-l. Re- 
spectively, in the N = 20 system with (Si, S 2 ) = (1, 4) 
the magnetization per cell changes from M = 0.5 to 
M = 0.4 (Fig. §). 




FIG. 3. Order parameters in the F and C phases as ob- 
tained from the SWT (thick lines) and the ED method for 
N = 20 with periodic boundary conditions (open circles). The 
thin broken line shows the ED data for the magnetization per 
cell Mo (N = 20). The classical values of the parameters are 
shown by dotted lines. 

More important changes take place in the C phase. 
The canted spin state is additionally characterized by 
the staggered transverse field (Sir) which breaks the 
spin rotation symmetry U(l) in the xy plane. In Fig. ^ 
we present result for the transverse order parameter as 
obtained from the SWT and the ED method. The ED 
data are obtained from the ED results for (Si*), by us- 
ing the relation (Si*) = (Sip) tan 0. The pitch angle 9 
is extracted from the spin-spin correlation function (see 
Fig. |)| 

(Sip • Slg 1+e 3 ) / q\ 

K Sl = -7= ~ p (9) 
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FIG. 4. The spin-spin correlators Ks 1 and Ks 2 vs. J2/J1 
used to obtain the pitch angle 6 between neighboring spins 
in the Si and 5*2 subsystems. The ED results (N = 20) are 
represented by open circles. The dotted lines show the re- 
spective classical values. The thick line represents the CCM 
results for cos 26 in the SUB4-4 approximation. 



which reduces to Ks t — cos 29 in the classical limit 
Si — > 00. The respective correlation function K$ 2 for 
the S2 subsystem can be defined in a similar way. Figs. ^ 
and H show that in the C phase the combined disorder- 
ing effect of the magnetic frustration and quantum spin 
fluctuations is strongly enhanced as compared to the F 
phase. SWT predicts a complete magnetic disordering 
of the S2 spins (m 2 = 0) starting at a = a*. The lat- 
ter-point exists in the SWT for arbitrary spins Si and 
SJ^ and preceeds the classical spin-flop transition point 
(a = 0.5). The effect is strongly pronounced in the ex- 
treme quantum case (1, i) where both the SWT and the 
ED data predict a* w 0.46. The extrapolated CCM data 
seem to show the same tendency. The steps in the ED 
data are due to finite-size effects and reflect the change in 
the total spin S t (a good quantum number) of the abso- 
lute ground state with the frustration parameter J2/J1. 
For instance, in the ED data a* appears as a point where 
S t changes from St = 1 to St = (respectively, M 
changes from 0.1 to 0). Although SWT formally predicts 
7711 > even beyond a* , the latter region is not accessible 
(for the standard SWT) as the phase with mi > and, 
7B2 = does not appear in the classical phase diagram.liil 
We suggest that the point a* is just the new position of 
the quantum spin- flop transition (see below). It is inter- 
esting to note that quantum fluctuations in the C phase 
do not change substantially the parameters Mq and 9 
(see Figs. | and §). 

Finally, the mixed-spin collinear state has properties 
which are similar to those of the collinear phase in the 
antiferromagnetic Ji — J 2 model, i.e., it is a Neel-type 
magnetic state characterized by the on-site magnetiza- 
tions mi,m 2 ^ 0. Both the SWT and the CCM predict 
that the on-site magnetization 777,2 vanishes at a point 
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(a c 2) which is very close to the classical spin- flop transi- 
tion point. 
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FIG. 5. Spin-spin correlators vs. J2 / Ji in the system with 
iV = 20 spins (ED results). The presented functions are scaled 
by their absolute values at J2 = 0. 

The ED data for different spin-spin correlations 
(Fig. ||) give some additional evidence in favor of the sug- 
gested phase diagram. Here the phase transition point 
a C 2 can be indicated as a point where the mixed-spin 
correlations (being ferromagnetic in the N phase) change 
their sign. We observe that close to the point a c2 the 
spin-spin correlations among the S2 spins are small and 
also change their sign. On the other hand, the corre- 
lations among the Si spins remain relatively large (and 
antiferromagnetic, as in the N phase) down to the point 
a* w 0.46, where the singlet ground state disappears. 
Thus, there are some indications that the established 
quantum spin phase (a* < a < a c2 ) has the symmetries 
of the canonical Neel state composed entirely of S\ spins 
(i.e., mi ^ 0, m 2 = 0). However, due to the finite-size 
effects in the ED data we can not exclude the spin-liquid 
state (mi = 0, r«2 = 0) as a possible ground state in 
the above region. Further understanding of the proposed 
phase diagram may be achieved by studying the low- lying 
excitations in these magnetic phases. 

C. Low-energy excitations 

In the F phase i?ik = -B 2 k = D k = and the disper- 
sion relations in Eq. ([5]) are simplified to 

45 = V^k-Ck± ~(,4i k - A 2k ) , (10) 

where A lk = 4JiS' 2 - 4J 2 Si(l - 1^), A 2k = 4Ji.5i - 
4J 2 52(1 - Vk), Ck = 4Ji V / S'i527k, v k = cosfcicosfc 2 , 
and 7k = (cosfci -I- cosfc 2 )/2. The low-energy physics 

is controlled by the acoustic branch u>[ F '. For small |k| 

u>\ takes the form 



where pf^ = J1S1S2 — ^($1 + Sf) is tne spin-stiffness 
constant of the ferrimagnetic state and mo = {S1—S2) /a, 2 
(a = v2) is the magnetization density. p s remains fi- 
nite at the F-C transition point a c \. The F-C phase 
transition is connected with a softening of the acoustic 
branch u>[ F ' at the corners of the PBZ. For instance, close 
to k = (0, 7r) the magnon spectrum reads 



J F > = 8Ji5i(a c i - a) + 4J x 8x [k\ + (tt- fc 2 ) 2 ] . 



On the other hand, the optical branch bj\ remains sta- 
ble at a c i and does not play important role in the quan- 
tum F-C phase transition. 
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FIG. 6. Spin-wave spectrum in the C phase along the e2 
direction of the PBZ (fc a = 0, Ji = 1). 

The excitation spectrum in the C phase is described 
by the general expression (H), where Ai k — AJ^S^t + 
4J 2 Si[l-i 2 (2^ k )], A 2k = 4J 1 S 1 t-4J 2 S 2 (l-z/ k ), B llt = 
-4J 2 Si(l - i> k , B 2k = 0, C k = 2JiV^(l + t)7k, 
D k = -2J 1 ^j%S 2 ~{l - i)7k, and t' 1 = 2aa (see Fig. §). 
As may be expected, in addition to the quadratic spin- 
wave excitations, Eq. (JTTJ) , close to the zone corners 
there appear linear Goldstone excitations connected to 
the spontaneously broken £7(1) symmetry in the C phase. 
For example, in the vicinity of k = (0, tt) the spectrum 
takes the form 

u[ c) = v^fkl + (Tr-k 2 ) 2 , v = AJ 1 S 1 ^a 2 -a 2 cl . 

Note that the spin-wave velocity v vanishes at the F-C 
transition point a c \. 
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FIG. 7. Spin-wave spectrum in the N phase along the e2 
direction of the PBZ (ki = 0, Ji = 1). 

Finally, the spin-wave spectrum in the mixed-spin 
collinear phase is given by the general relation (||), 
where A lk = 4J 2 5i, A 2k = 4J 2 5 2 , B lk = 4J 2 Si I / k , 
-B 2k = 4J 2 5' 2 i/ k , C k = 2Ji^SiS 2 cos fc 2 , and L> k = 
2JiVS'i5 2 cos k\ (see Fig. 0). As compared to the an- 
tiferromagnctic J\ — J 2 model, the excitation spectrum 
in the mixed-spin collinear phase shows some important 
peculiarities which specify the instability at the classical 
transition point a = 0.5. First, in the mixed-spin sys- 
tem the degeneracy of the spin-wave branches aj[ N ^ and 
ct>2^ is completely removed. Second, the spin-wave exci- 
tations described by u>2 formally remain stable down to 
the point a = y/a/(a + 1) [= \/2/3 « 0.47 for the (1, ±) 
system] . The phase transition at a c2 is entirely driven by 

instabilities in the lower oj^ excitation branch. At the 
classical spin- flop transition point (a = 0.5), the branch 

u>x exhibits two full lines of zero modes (i.e., k± = 
and fc 2 = 0) in the PBZ. The soft lines are connected 
with ferromagnetic fluctuations in the ei and e 2 direc- 
tions on the square lattice, so that we can expect that 
quantum fluctuations do not change the position of the 
classical transition, point (in accord with the data pre- 
sented in Fig. §)]i3 



III. CONCLUSIONS 

In conclusion, we have studied the quantum spin 
phases in the mixed-spin Ji — J 2 Heisenberg model de- 
fined on a square lattice. The reported results sup- 
port a phase diagram with an additional quantum spin 
phase close to the classical spin-flop transition. In par- 
ticular, in the extreme quantum case (1, g), we have 
found indications of a quantum spin phase (in the re- 
gion 0.46 < a < 0.5) which is characterized by magneti- 
cally disordered S2 spins and Neel long-range ordered Si 
spins. The latter suggestion is based (i) on the analysis 
of the low-lying excitations showing that the mixed-spin 
collinear state survives at most down to the classical tran- 
sition point a = 0.5 and (ii) on the SWT and ED results 



for the sublattice magnetization in the C phase m 2 : both 
methods predict that m 2 vanishes at a point (a ~ 0.46) 
preceeding the classical transition point a = 0.5. How- 
ever, due to strong finite-size effects in the ED data and 
the qualitative character of the SWT analysis, we can not 
exclude the spin-liquid state as a possible ground state in 
the above region. Studies with other analytical and nu- 
merical methods would be necessary to characterize the 
true ground state in the above region. 
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